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Abstract —In this paper, a new model for traffic on roads 
with multiple lanes is developed, where the vehicles do not 
adhere to a lane discipline. Assuming identical vehicles, the 
dynamics is split along two independent directions—the K-axis 
representing the direction of motion and the A-axis representing 
the lateral or the direction perpendicular to the direction 
of motion. Different influence graphs are used to model the 
interaction between the vehicles in these two directions. The 
instantaneous accelerations of each car, in both X and Y direc¬ 
tions, are functions of the measurements from the neighbouring 
cars according to these influence graphs. The stability and 
equilibrium spacings of the car formation is analyzed for usual 
traffic situations such as steady flow, obstacles, lane changing 
and rogue drivers arbitrarily changing positions inside the 
formation. Conditions are derived under which the formation 
maintains stability and the desired intercar spacing for each of 
these traffic events. Simulations for some of these scenarios are 
included. 

Index Terms —lane-less traffic model, formation control, 
multi-agent system 


I. Introduction 

This paper is motivated by the desire to develop models 
to analyze traffic when the roads are wide and the lanes 
are blurred. Such traffic behavior is not uncommon in many 
roads in India where, to misquote J K Galbraith, a ‘func¬ 
tioning anarchy’ prevails. We develop a stylized model to 
describe the local interaction among the vehicles on such 
roads and use this microscopic description to characterize 
the emerging macroscopic behaviour. In effect, we consider 
a multilane system except that we assume that there is no 
strict demarcation of lanes and that vehicles are affected by 
those in a ‘cone’ rather than by the vehicle right ahead. On 
such roads, multiple cars drive abreast but do not adhere to 
lane discipline. Toward this we adapt the single lane model 
of [6] to develop our multilane model. Specifically, a directed 
graph is used to model the influence on the acceleration of a 
vehicle by others in its neighbourhood, i.e., those that are in 
its cone. We seek an equilibrium analysis of the dynamical 
system model that we develop. Our notion of stability refers 
to the condition that all cars attain the same velocity as the 
leader. Since this is a car following model, we also have the 
notion of levels and the analysis will also obtain conditions 
such that cars in one ‘level’ maintain a fixed spacing from 
cars in the next ‘level.’ Our analysis will primarily use the 
Laplacian of the directed graph that models the influences, 
which in turn will allow us to dissolve the lanes, so to say. 


The influence graph can also be a weighted graph to enable 
us to model the relative degrees of the influences. 

Models for vehicular systems have been widely studied 
both from a microscopic [5] [7] as well as from a macro¬ 
scopic [9] [22] perspective. Various forms of stability such 
as input-to-state [30], string stability [29] and mesh stability 
have been considered. Mesh stability ensures damping of 
perturbations in vehicle formations, as they move away from 
the source [27][21][28]. Various applications are studied 
namely obstacle avoidance in [20] and lane changing in [14] 
[8]. In [11] multilane model for vehicular traffic is considered 
and a fluid model is developed. Unlike string or mesh 
stability theories, where the main concern is magnification 
of disturbances down infinite chains of cars, in this paper, 
we study finite number of cars. While it is relatively easy 
to show stability properties in this case, we are concerned 
primarily about equilibrium spacing between the cars. We 
show that it is possible, with purely distributed control laws, 
to achieve and maintain desired spacing, even in an ad-hoc 
laneless traffic situation under a variety of common traffic 
conditions and disturbances. 

The results presented here, as well as the tools used to 
derive those results, are heavily dependent on the theory of 
consensus in multiagent systems (see e.g [4] [10] [19] [24] 
and the references therein). While [26] has studied vehicle 
consensus, [23] [17] derived stability conditions for time 
varying topologies. Most of our analysis depend on these 
recently developed theories and particularly on [25], where 
arbitrary vehicle formations (not necessarily unidirectional as 
in roadways) have been studied. We specialize these results 
for ad-hoc laneless road traffic, and in the process derive 
distributed control laws that preserve inter-vehicle spacing 
and stability under time-varying formations, switched in¬ 
fluence graphs and impulse changes in driving objectives. 
While these properties are hard to derive for general types 
of motion, it becomes relatively simple here due to typical 
uni/bi-directional constrained motion possible on roadways. 

The rest of the paper is organized as follows. In Section [II] 
we set up the notation and our assumptions on the models 
for the X and the Y directions. We then describe the control 
laws that are used by the vehicles for each of X and Y 
directions. In Section[IIIJ we analyze the dynamics of motion 
in the Y direction and obtain a stability criterion. We also 
characterize the influence graph. A similar analysis is carried 
out for the X direction in Section IV In Section[V]we present 



the analysis when the influence graph varies with time. In 
Section VI we analyze the effects of impulses on the stability 
of the formation. Finally, in Section |VII| we present some 
simulation results. 


II. Traffic Model and Control Laws 

As discussed in the introduction, we are interested in 
analyzing the motion of a finite number of cars along 
roads which has no lane demarcations. In other words, cars 
can, in principle, occupy arbitrary positions with respect to 
other cars as long as they do not collide with each other 
or with the side of the road. However, all the cars want 
to travel along the road in one direction and would also 
like to keep safe a distance from all its neighbours based 
on visual feedback about the position and velocity of the 
neighbouring cars. We name the direction of travel along 
the road as the Y -axis and the direction perpendicular to 
the road as the X-axis. Though in reality, a car’s ability to 
maneuver in these two directions is coupled, for simplicity 
we assume that the dynamics of each car along these two 
directions are independent. Moreover, the control laws for 
each of these directions consider a different set of influencing 
neighbours. This later assumption is realistic since a driver 
typically looks to the side before moving sideways, while 
he considers only the cars roughly in front of him for 
normal forward motion. In effect, the influence graphs for 
the X and Y directions are different in our calculations. 
In addition, we also assume, again for simplicity, that the 
cars in the formation are identical. We state the various 
modeling assumptions along these two directions after a brief 
description of graph theoretic notation. 

A. Notation 

The sets of naturals, reals, positive reals and real n tuples 
are denoted by N, R, R + and R”. An undirected graph 

= ( Y,S,w ) is a finite set of nodes f connected by a 
set of edges § C V x "Y along with a function w : § —> R + . 
When two nodes a, £ 'Y and aj £ 'Y are connected to each 
other the graph said to have an edge between a, and cij , 
denoted by (a,,a ; ) £ S. A graph is said to be connected 
when there exists a path between any two nodes. A spanning 
tree is a connected graph %\ ree = {'Ytree,<%tree) having no 
cycles in the graph. If the edges of a graph 'Y = {'Y. Y) 
are directed i.e. (a,-,a,-) £ $ =£> (a,-,a;) £ S, the graph is 
called a directed graph (or digraph). A rooted directed tree 
is a digraph such that there exists a node (called root) and 
a directed path from that node to all other nodes in the 
digraph. A digraph r Y is said to contain a directed spanning 
tree, if there exists a rooted directed tree ree — {Y tree ,<f t ree) 
such that 'Yin,,, = "Y and S tree C S. We define the Laplacian 
(L) for a directed graph with weights w,y as follows: := 

— 1 if and £u = YJj=\ w ij := indegree. 

B. Assumptions for Y -direction 

Al) For determining the influence in Y direction we 
consider each driver to have a fixed viewing angle ( aov) of 
120° and a car must be within this viewing angle to influence 



(a) Y direction motion: car 3 is (b) X direction motion: black lines 
not seen by car 1 indicate the underlying graph 


Fig. 1: Different influence graphs for X and Y direction 
depending on the angle of view (aov) of the car 


the driver. This forms the conical section c/itLcL as shown 
in Figure [Ta] In this figure cars 2 and 4 lie within the conical 
section d^dj formed by car 1 and hence can influence car 
1. Car 3 being outside the cone cannot influence car 1. 

A2) We assume that the convoy follows a fictitious leader 
in the direction of travel. As discussed later, the role of the 
leader node is to set a desired velocity for our formation. 
This fictitious leader can be considered as a mathematical 
representation of velocity limiting rules, regulations or driver 
experience in practice. In Figure 3a node 0 is the leader. 

A3) The influence graph for Y axis is assumed to have 
a directed spanning tree with node 0 as the root. This 
assumption is equivalent to connectedness of the influence 
graph in this framework. 

A4) Cars with same number of hops from the leader node 
are referred to as being in the same ‘level’. In Figure [3a| cars 
1,2, and 3 belong to level one. 

A5) We assume that influence of any car extends up to 
a maximum of one level for ease of exposition. This is 
elaborated on in later sections. 

A6) The convention used for numbering cars in the graph 
is as follows: (i) Cars are numbered according to their levels 
with cars in higher levels having higher numbers, (ii) Cars 
in the same level are numbered from left to right with left 
most car being the highest number in that level. 

A7) We assume that each node can choose the distribution 
of weights on its incoming links. This is related to the 
driver’s discretion about relative importance of the various 
cars in front of him. 


C. Distributed Control Laws for Y-direction 

Consider a convoy of cars represented by influence graphs 
as shown in Figure [3a] along Y direction, where node 0 is 
assumed to move at an arbitrary velocity v'q. We assume that 
in the direction of travel, the last car has the lowest nonneg¬ 
ative y coordinate: car 9 will have the lowest nonnegative y 
coordinate and car 0 has the highest y coordinate in Figure 
[3a| Let the velocity and position of car i for motion in In¬ 
direction be denoted by vw and y t . The proposed control law 
model for car i in Y direction is a modified version of single 
lane car following model of [6] and is as follows, 

yi = Vyi 

Vyi’-= Yj b ( w ij v yj- w ii v yi)+ k (y v iiyj-Wiiy i -gy) ( 1 ) 






Here jV (i) represents the neighbour set of car i, w,j is 
the weight given to link connecting node j — » i V j ^ i, 
wu = Yj"j= i w ij —■ W, b and k are constants, and parameter 
g y is a tuning parameter used to adjust the equilibrium 
distance between cars in consecutive levels. Consider an n 
car formation with node 0 as leader moving along the road. 
Let y £ R' 1+1 represent the y coordinates or absolute position 
of n + 1 cars. v y €K" +1 represent the velocities of the cars. 

A graph = ('V, S, w) having a node set V and an edge set 
S' is used to denote the influence diagram between various 
cars as per Al. Each node in the graph represents a vehicle 
(agent). The directed edges are introduced as follows: If 
vehicle j can be sensed by vehicle i (i f j), then an edge 
from j to i exists with weight w,y and is denoted by j —> i. 
Then the control law in 0 can be written as follows: 
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( 2 ) 


where leR”, 0 € M" , V denotes the Laplacian of the 
directed weighted influence graph. 

D. Assumptions for X-direction 

Bl) For motion in A-direction the angle of viewing 
extends to 180°. This is reasonable since a driver would 
usually look sideways before making a sideways deviation. 
In Figure [Tb] car 1 is influenced by cars 2,3, and 4. 

B2) For horizontal motion, the boundary of the road 
needs to be incorporated into the model. This is achieved 
by assuming some fictitious cars are moving along the 
boundary of the road at each ‘level’. Levels are still defined 
using the T-direction influence graph and hops from the 
pseudo 0-node. In Figure [3b] cars 1,5,9, enclosed by the red 
block, represent the boundary of the road. These nodes have 
same dynamics as the regular cars but move under slightly 
specialized laws: (i) these nodes have no horizontal velocity 

(ii) they continue to have the same y-coordinate and the 
same y velocity as the leftmost car in the corresponding level 

(iii) boundary nodes are influenced only by other boundary 


nodes from the levels directly above them. In Figure 3b 
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Fig. 2: The convoy of cars for which the influence graphs 
are represented below 




(a) Information flow along 
Y axis 


(b) Information flow along 
X axis 


Fig. 3: Figures show the influence graphs for both directions 


where ,yL(() represents the neighbour set of car i, Wij is 
the weight given to link connecting node j —> i V j 7 k i, 
wu = Y!}= 1 w ij =: W, b x and k x are constants, and parameter 
g x is tuned for desired equilibrium distance between cars 
in the same level. This law differs from the T-axis law 
in that it requires some extra parameters z,’s to ensure the 
desired inter-car spacing. The X axis motion of n cars in 
the formation is described below. Let x £ R" denotes the 
coordinates and v x £ R" represent the velocities in the X 
direction. The control law in 0 is as follows: 


(4) 
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boundaries on both sides can be incorporated in the model, 
here we analyze the system with boundary on one side of 
the road for simplicity of exposition. 

B3) For X axis motion node 1 is assumed to play the role 
of the leader node as shown in Figure [3b] This again sets 
the x velocity of the formation. 

B4) The influence graph for X direction is assumed to 
have a directed spanning tree with node 1 as root. 

E. Distributed Control Laws for X-direction 

Let the velocity and position for car i be denoted by v X i 
and Xj. The proposed control law for car i in A-direction is 

Xi = v xi 

v X i := b x (wijV xj -wuv x ^) +k x (wijXj - Wuxi - Zig x ) 

J=-A' ( i) 

(3) 


nxn , c £: 

0 . 


where leR" and 0 £ R" and 0 £ 

C:= [0 Z 2 ■■■ Zk 0 Zk +1 u ... z n \ 

where the nonzero z \, 7.2 - ■■ ■ z, n are constants locally com¬ 
puted by each car to ensure a spacing of g x between cars 
in the same level. The L x denotes the directed weighted 
Laplacian of the A-axis influence graph. The dependence 
of Zi and g x is elaborated on in later sections. 

III. Analysis of T-dynamics 
In this section we give conditions under which control law 
0 gives the desired spacing with unidirectional communi¬ 
cation for Y axis motion. We rewrite 0 with the position 
(yo(0) anc * velocity (vo(l)) of leader node 0 as external 
inputs: 
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where B y £ 


i>nx2 


y, W G 


L y is the reduced Laplacian 


links from car j\ to car i be denoted by w tJI 


obtained from L y after removing the row and column corre¬ 
sponding to node 0. B y contains the columns from L y which 
denote the links between the leader node 0 and the remaining 
cars in the formation. For Figure [3a] assuming unit weights, 
B y is given by. 


B y = 


0...0 

0...0 


1 T 


( 6 ) 


In (|6|, the l’s in B y represents information flow from node 
0 to cars 1,2,3. We give our first main result for Y direction 
motion of a convoy having node 0 as leader. Let 


y-=[y w] J 


and T y := 


0 I 

-k V - b V 


and so 

on. The i ,h row in ([8]) can be expressed as follows: 
Wijjj i + w ih y h ...+ w ijn y jn = g y + Wy t . From induction 
hypothesis, yj = yi+ g y /W. This completes the proof. □ 


Clearly, for unweighted influence graphs, the conditions 
for Theorem 0 reduces to each car having the same 
indegree (say M), thereby guaranteeing the relative spacing 
between two consecutive levels to be g y /M. 

IV. Analysis of A-dynamics 
In this section we analyze the vehicle motions in X- 


(7) 


Theorem 3.1: Consider a weighted directed graph such 
that the total weight ( W ) across all incoming edges is the 
same for each node. Then the autonomous system y = T-'y 
achieves an asymptotically stable equilibrium point at the 
origin (0 £ R 2 "). Moreover if the leader velocity v y o is 
constant, then: 

1) Vyi -> v v o as t °° V i e {1,.. .«} 

2) \yi(t) — yj(t ) | — > 0 as t -A °° for all i, j in the same level. 

3) At equilibrium, the relative spacing between cars in 
two consecutive levels is g y /W. 


direction under the assumptions stated in Section II-D The 
analysis, while being similar to the Y -dynamics, differs in 
some crucial features. One of the main differences arise since 


the influence graph can be bidirectional (see Figure 3b i: in 


other words, two drivers can simultaneously look towards 
each other and decide on their A-control. Secondly, recall 
that the A-velocities of the pseudo (road boundary) cars 
(e.g. cars 1, 5, 9) are zero by assumption (e.g. v x \ = v A -s = 
v .\9 = 0), and hence their X-positions are always constant 
(xi = xj = xg = constant). Hence when we rewrite 0 with 
node 1 velocity and position as external input, 
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leader, has non-negative real eigenvalues and 0 is a sim¬ 
ple eigenvalue, with [l 0 ... 0] as its left eigenvector, 

and the Laplacian right nullspace consists of the vector 
[1 1 ... 1 ] 7 . Additionally the following result holds. 

Lemma 3.2: Under the numbering scheme in Assumption 


where x, v x £ R" , B x contains the columns from L x which 
denote the links between the leader node 1 and the remaining 
cars in the formation. We assume that every car wants to 
position themselves at a distance g x from both adjoining cars 
in the same ’level’. Let, 


II-B the influence graph has no cycles and the Laplacian 


( U ) has a lower triangular structure. Moreover, the diagonal 
entries of V will be same for all rows excluding the rows 
representing level one and node 0. 

Proof of Theorem |J. /[ 

We first analyze the asymptotic stability of the autonomous 
part of the system in 0 From [25] we get Re(kfT)) < 
0 V i, V b, k > 0, if A/(-L- v ) < 0 V j and lj(L y ) £ R V j, 
where Re (.) denotes the real part. The equilibrium point 
for the autonomous part of the system in ((5) is v y = 0 and 
- L y y = gy 

Next we analyze the system in ([5]) at equilibrium with 
fixed velocity vo of node 0 as input. The required spacing 
can be computed from. 


x:=\x v_y] 1 


and T x := 
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-k x L x 
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-br L X 


( 10 ) 


Theorem 4.1: Under the above assumptions, there exists 
h, k such that the autonomous system x = F‘x asymptotically 
achieves equilibrium at the origin (0 £ R 2 "). Moreover, 

1) Vxi(t) — > 0 as t — > °° V i e {1,... ,n}. 

2) There exist C such that |x; — x,-+i| —> g x as t —> °° for 
i,i+ 1 in same level. 

3) The C achieving (2) can be computed locally. 

For obtaining a spacing of g x between cars in the same 
level, we impose the additional constraints for cars i and i +1 
in the same level (denoted by J'') as follows: 


-L y y = g y 


( 8 ) 


Xi+ 1 Xi — g x Vi, (+16/ 


( 11 ) 


where 1 £ R" for n + 1 nodes. Now we use induction. Cars 
in the first level will have one link connecting to the node 
0 by A2. For any car i in the first level from ([S} we get, 
yo = yi + gy/W . This will ensure that all cars in the first 
level have the same y coordinate. 

Consider some car i in level . Let cars in the level 
above be denoted by j\. jj, ■ ■ ■ -jn- Let y coordinates 
of the cars be same. Let the weights on the 


Proof of Theorem \4.1\ The stability condition on constants 
b and k in control law 0 can be easily obtained from [25]. 
Claim (1) follows immediately. 

Next we show the existence of C = [0 zi ... z„] T in 0 
for cars in level JT. The same arguments can be repeated for 
other levels. Let m and m be the total number of cars in level 
y and in the level above {J* — 1). Let the cars in level J/ be 
linked to cars jk with weights given by Wjj k V k £ (JT — 1). 

























Then the following equation can be obtained from Q at 
equilibrium: 

L x x = -g x C (12) 

Choosing rows corresponding to cars in level from Cl- 
and combining with (fTTj), we get 


" 1 

0 

-1 0 


0 

... 0 " 


f Xj -j 


r o i 

0 

w Ui 

wii wu+i 

-1 

0 

... 0 





0 

w MJi 

... wi+li+l Wi+u+2 

0 

-1 

... 0 


x m+m 


6 








gxZ\ 


gx 

0 


-1 1 

0 

0 

0 0 





0 


0 -1 1 

0 

0 

0 0 












. 8xZm+m . 


.gx. 


(13) 


We show the existence and uniqueness of constants Zi using 
rank conditions. Rearranging columns of the augmented 
matrix from © we get: 
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The matrix in ( fl4| ) is in row echelon form with unity pivot 
elements. Thus (jT4j is full row rank matrix. This guarantees 
the existence of at least one solution ( n ) giving the desired 
equilibrium spacing. 

The local/distributed nature of the C computation 
can be verified by considering © for car i: 
wtjiXji + ... + w ijn Xj„ - waXi = gxZi and x l+ \ - x t = g x . 
Clearly these equations only require information from the 
immediate x-neighbours of the car i. □ 


Remark 4.2: At the expense of some new notation, an 
explicit formula for C can be given easily. Let xr £ R" be 
such that Xfg x £ R" give the x-coordinates of all the cars 
in the formation. For example in Figure 3b if we want 
all cars in the same level to have maintain a spacing of 
g x from each other the corresponding xj should be, xj = 
[0 1 2 3 ...0 1 2 3} T . Then it is easy to verify that the 
constant vector C in control law {4]i is given by C = Lxf. 

Remark 4.3: Though the above results have assumed that 
influences percolate across atmost one level, this assumption 
can be easily extended to multi-level influence graphs. Such 
a requirement might arise, e.g. for a driver looking in the rear 
view mirror while taking turns to avoid cars behind him. 


V. Time varying graphs 

In this section we analyze situations where the influence 
graph changes over time, due to relative motion of the 
vehicles during various traffic events, e.g., road widening or 
traffic signals. During road widening scenario, cars moving 
in a convoy may spread out. As the cars approach traffic 
signals the cars in the formation may come closer. As this 
motion takes place, cars may move in and out of influence 
cones of other agents, thereby changing the influence graphs. 


the transition matrix and the Laplacian. Hence {2]) and 0 
become switched systems. We assume, quite reasonably, that 
this switchings occur after finite intervals of time. Let P 
be a finite index set: P = {1,2Let <7 be piece wise 
constant switching signal <7 : [0, °°) P, r<j be the transition 
matrix and L a the graph of the Laplacian of the system in ([2]) 
and ([4]) corresponding to cr(f). We analyze the Y -dynamics 
only for lack of space. Similar arguments work for the X- 
dynamics as well. The switched version of 0 is given by. 
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(15) 

where B y a £ M” x2 , y.v y £ R". Clearly, L y a . L y a , and B y a 
changes with the switching influence graph. In this section 
we assume that the cars still want to preserve the same inter¬ 
vehicle spacing even though the graph is changing. Hence 
the spacing constant g y remains same for all (7(f). We show 
that this is possible by local re-computation of the weights 
assigned to the incoming edges by each car. We now give 
the main result. 

Theorem 5.1: Under the assumptions stated in Section II- 
|B| the autonomous part of system ( p~5j ) is globally uniformly 
asymptotically stable and © is BIBO stable. Moreover, if 
v-yo is constant and the net weight of all incoming edges for 
each node be kept same (IV) for all <7(f), then: 

1) v yi — > Vyo as t — t °o V i £ {1, ... n} 

2) |>’;'(f) — yj(t ) | — > 0 as t —> °o for all i, j in the same level. 

3) At equilibrium, the relative spacing between cars in 
two consecutive levels is g y /W. 

Proof of Theorem |5. 1 [ [Sketch of Proof] Stability for each 
subsystem in ( p~5] > was shown in Theorem 3.1 The transition 
matrix for the system subject to switching signal cr(f) in © 
is given by 


r y — 

1 a — 


0 7 

-kL y a —bL y a 
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where R = R T > 0, I = I n xn/<T q is a constant satisfy¬ 


ing l/kq < Amin (La). Noting that L y J + L y a > 0 for q as 


stated here 
Section II-B 


it is eas 


isy to verify under the assumptions of 
that rjj P + PT y a < 0. It follows [15] that the 
autonomous part of ( p~5j ) is globally uniform asymptotically 
stable (GUAS). Further it is well known [16] that for a linear 
switched system, GUAS ensures BIBO stability. Claims (1), 
(2) and (3) follow in a similar fashion as in Theorem 1 3. 1| 
Remark 5.2: We have assumed that the net weight for 
each node remains same. At every switching instant, we 
assume that each car can locally redistribute the available 
weight on its links. With this the equilibrium point of 
the autonomous system in © remains invariant across 
switching. 


A. Obstacles 

Using the switching theory developed above, stability of 
inter vehicle spacing in the presence of stationary obstacles 
can be ensured. The obstacle is modeled as a stationary 
car having zero velocity in both directions and is numbered 
































(a) Graph before 
obstacle is sensed 



(b) Graph when 
obstacle is sensed 



0 


(c) Graph afer 
obstacle has passed 


Fig. 4: Changes in information flow graph for Y direction 
when 4 cars encounter an obstacle denoted by node 1 


consecutively along with other nodes. Figure [4] shows the 
transitions in the influence graph for a typical obstacle 
crossing. We briefly describe the dynamics in Y -direction. 
Similar arguments hold for the X-dynamics. 

Let i’o. v, and yo. y, be the velocity and position inputs of 
leader node 0 and obstacle denoted by node i respectively. 
By hypothesis, v,- = 0 and y, = constant. The control law for 
motion in Y direction: 



0 

-kL y c 


I 

y 


0 


0 

bL y a 

. v y. 

-kgy 

1 

+ 



bv o 
kyo 
bvj 
kyi 


(16) 


where B a £ M"~ lx4 , y, v v £ R' 1-1 . As usual, B a identifies the 
(switching) links between leader node 0 and the remaining 
cars (except node i) and between node i and the remaining 
cars in the formation. 

Corollary 5.3: Suppose all cars influenced by the obstacle 
are also influenced simultaneously by at least one more car 
in the formation. Then Theorem |5. 1 1 holds for (fl6|). 


B. Lane Change 

In lane changing scenarios we assume that a particular 
agent decides to change his X position arbitrarily within 
the formation. We assume that the maneuver takes finite 
time to complete. This is illustrated in Figure [5] where car 
6 is changing lanes within the formation. Suppose car i is 
changing lanes. We assume car i to be an external input. The 
Y -axis dynamics is given by: 
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(17) 


where y,y, v y ,v y £ M" -1 , y = y + g y . As long as a directed 
spanning tree continue to exist during the lane change, 
clearly Theorem 5.1 continues to hold for O- Similar 
analysis holds for the X dynamics but is not included for 
lack of space. 


VI. Stability under impulse effects 
In previous sections the individual cars were kept equidis¬ 
tant from each other using the constants g y , g x , xj (and 
C). However the cars might want to change the equilibrium 
spacing in a variety of situations, such as, when the road 
narrows/widens and we want cars to come closer/further 


or some rogue drivers want to change their position in 
the formation arbitrarily. Individual cars can decide on a 
new g y , xy or C suddenly, introducing discontinuities in ([2]) 
and ©• Such discontinuities are known as impulse effects 
[12], For simplicity we assume that the influence graph 
remains invariant and the velocity vectors are continuous for 
t £ [0,°°). However, time varying graphs can also be handled 
easily. 

The autonomous part of the system in 0 with impulse 
effects at fixed time instants Tj. (assume 0 = To < T] <Xi < 
... < lim^oo T/; = °o) can be characterised as follows, 
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where y = y + g y . Let the cars decide to change instante- 
neously the spacing constant from g y (t) to g y (t + ) = g y (t ) +g y 
at t = Xk- This creates a discontinuous change in the position 
state vector y(f + ) = y{t) + Ay where Ay = gy at t = ro¬ 
under these assumptions, the existence and uniqueness of 
solutions for the system in ( fT8j ) is ensured [32, Theorem 
1 . 6 . 2 ], 

Theorem 6.1: Suppose L remains unchanged for t £ [0,°°) 
and Vy(f + ) = v y (t). Then ( fl8| is exponentially stable if g y (t + ) 
is chosen such that llvff + )ii<ibm 


Proof of Theorem 6.1 Let y := [y v y ] . The common 


Lyapunov function at time t Xk is given by, V(t,y ) = 


I/q 

0 


y where 4^ < A mm (L). At the switching instant 
t = x k , V(t + ,y + g y ) = (;y + g y ) T !. ® (y + lv) where g y := 


\g y 0 ] , gy £ ]R 2 ". The 0 £ g y exists because we have 
assumed that the velocity vector is a continuous for t £ [0,°°). 

From [12], and [32, Theorem 4.14]: V(t,y) > V(t + ,y+g y ) 
guarantees exponential stability. It is easy to be verified 
that this condition is satisfied if ||y(f)|| > ||y(f + )|| This 
completes the proof. □ 


A similar result holds for impulse effects in X -dynamics. 
VII. Simulations 

In this section we present some numerical simulations. In 
our numerical simulations we use the scenarios of changing 
formations, obstacle avoidance, and lane change. We use the 
initial conditions depicted by graphs in Figure [6] for Y and X 
axis motion respectively. Unidirectional communication for 
X and Y axis is assumed with an influence level of one. 
Cars 1,5,9,13 indicate the boundary of the road in Figure 
[6] The maximum number of cars in each level is 4. The 
angle of view for Y axis is a 120°cone and for X axis it is a 
180°angle. The system parameters in ([2]) and are set as 
follows: b = 0.4 and k = 0.001. The vertical spacing constant 
between consecutive levels is g y = 50 and horizontal spacing 
for cars in the same level is g x = 30. 

A. Changing formations 

When the equilibrium point of X or Y dynamics is 
changed we observe a change in the formation. We change 
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Fig. 7: The X vs Y coordinates of cars are indicated at a 
time instant just before changing formations. The system is 
at equilibrium at this time 
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Fig. 8: The X vs Y coordinates of cars are indicated at a 
time instant after formation has changed and the system has 
settled at equilibrium 
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Fig. 9: Figure shows the formation at the instant before car 
6 changes lanes 
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Fig. 10: Figure shows the formation at the instant after car 
6 has changed lanes 


vector Xf in Remark 4.2 at different time instants to attain 


differing spacing between cars in the same level. Changing 
Xf introduces impulse effects in the system in ([2ji. Figures [7] 
and [8] demonstrates this effect. Cars 1 to 4 act as references 
and are spaced 30 units from each other. Nodes 1,5,9,13 
represent the boundary. The V and ‘A’ represents the system 
a snapshot of the system at time 2000 and 5000 units 
respectively. The corresponding values of xj for which the 
figures are obtained is as follows. For time units of 5000 
units, x T f = [ 0 1 2 3 0 0.5 1 1.5 0 1 2 3 0 0.5 1 1.5]. For time 
units of 2000 units, xj = [ 0 1 2 3 0 1 2 3 0 0.5 1 1.50 1 23]. 

B. Obstacle avoidance 

Here we show the impact of a stationary obstacle on the 
convoy. In our simulation the obstacle affects all cars below 
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Fig. 11: Trajectories of cars in X direction changing due to 
the obstacle 
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Fig. 12: Velocities in Y direction are plotted as function of 
time. The deviation is due to the obstacle 

cars in level 2, i.e. cars 10 onwards. Figure [TT|plots deviation 
of x coordinates with respect to time. In Figure |TT] the red 
star denotes the stationary obstacle. Car 12 goes around the 
obstacle. This causes a change in position and velocity of 
cars 11 and 13 in the same level. Note that car 11 has moved 
towards the obstacle, this is caused due to the control law 
forcing car 11 to maintain a fixed distance with respect car 
12 and the obstacle. 

Figure [12] shows the change in Y velocities due to the ob¬ 
stacle.The black dotted trajectory is the reference trajectory 
in the absence of an obstacle. As cars 11 to 13 sense the 
obstacle their velocities drop as shown in the Figure. This 
drop causes a deviation in the y coordinates of the cars. We 
now present the case for a particular car arbitrarily changing 
x coordinates within the formation. 


C. Lane change 

For this example, car 6 changes lanes from the left to the 
right in the formation. In doing so car 6 has moved from 
one level below to one level above in the formation. Figure 
[9] shows the formation before lane change. As car 6 changes 
lane other cars are also affected. This can be seen in Figure 


10 in which cars 7,8 are moving closer to car 5. Cars 10,11 


are moving closer to cars 5,7,8 in Y direction. Note that xj 
has not been changed for level 2 even after car 6 has moved 
out hence the spacing is remaining same as before. 
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